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In the recent series of papers (cond-mat/0402471, cond-mat/0403618, cond- 
mat/0407618, cond-mat/0501586), Janis and Kolorenc discussed the role of 
the diffision poles in the Anderson transition theory. Their picture contra- 
dicts the general principles and is shown below to be completely misleading. 
Correct setting of the problem is given and the contemporary situation is 
discussed. The critical remarks are given on the relation of the diffusion 
coefficient with multifractality of the wave functions. 

The present paper has an aim to discuss the controversy in the literature concerning 
the role of the diffision poles in the Anderson transition theory. 

1. The spatial dispersion of the diffusion coefficient D(u,q) has a fundamental signifi- 
cance for the whole Anderson transition theory. 
It is well known (see e.g. [1, 2]) that a quantity 

0^(r-r') = (G| + Jr,r')G^(r'r)) (1) 

in the momentum representation has a diffusion pole 

= ■ 2 ™nf ] S 2 + treM (2) 
—iu + D(ui,q)q 2 

where D(uj,q) is the observable diffusion coefficient, G R and G A are nonaveraged retarded 
and advanced Green functions for an electron in the random potential, v(E) is a density 
of states for the energy E. In the localized phase, the diffusion pole transforms into the 
Berezinskii-Gorkov singularity 

t RA ( q ) = ^^A(q) + <j )re9 (q), (3) 

A(q) = J dve-^A(v) , A(r) = T^y (E |^(r)| V s (0)| 2 5(£ - e s )"j , (4) 

where ip s {r) and e s are eigenfunctions and eigenenergies of the electron in the random 
potential. In the original paper [3], the singularity (occuring from the terms with 
s = s f ) was established for the density correlator 

S(r -v') = -^ (jT, r s (r)Mr'WAr')M*)KE - e s )5(E + u>- e s ^ (5) 
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but it can be easily transform into the 1/uj singularity for <p RA [4, 5, 2]. Comparison of (2) 
and (3) shows that D(u, q) ~ uj in the localized phase: a slower dependence would destroy 
the 1/uj singularity in Eq. 3 and a more rapid dependence would lead to disappearance of 
the g-dependence in the singular part of (2), which surely exists according to Eqs. 3,4. As 
a result, 

D(uj,q) — —iuj d(q) , (6) 
where the limit uj — > is taken in the function d(q). It is easy to see from the relation 



1 + d(q)q 2 



A{q) = J dve- tq r A(v) (7) 



that d(q) is a regular function of q 2 , since all coefficients in the expansion of (7) in the powers 
of q 2 are finite due to exponential decay of A(r) at large r. The analogous regularity in q 2 
is expected for D(0,q) in the metallic state, while the anomalous spatial dispersion of the 
type q a can occur at the transition point. 

One can see from Eq. 6, that D(0,q) = in the localized phase; so the Anderson 
transition does not reduce to vanishing of D(0, 0) but has essentially more deep nature [6]. 
The question arises, how the spatial dispersion of D changes near the Anderson transition. 
One can suggest, that D(0,q) vanishes at the transition point simultaneously for all q. 
Such possibility looks incredible from viewpoint of phenomenological considerations in the 
spirit of the Landau theory. Indeed, all expansion coefficients of D(0, q) in q 2 should vanish 
simultaneously, irrespective of the way the critical point is approached and its location 
on the critical surface. Obviously, such vanishing cannot occur by chance, and should be 
backed by some profound symmetry. Does such a symmetry exist, and what is its origin? 
If this symmetry is taken for granted, then the order parameter for the Anderson transition 
should have an infinite number of components. Another possibility suggests, that D(0,q) 
vanishes for a single value of q, and then instability develops (in the spirit of a soft mode) 
leading to a first-order phase transition. In this case, one has to suggest an appropriate 
scenario. One can see, that even a type of the Anderson transition cannot be understood 
without solving the problem of the spatial dispersion of the diffusion coefficient. 

The situation is aggravated to the utmost by the existence of the Ward identity [1] 

AS k(q) = ^E^ kl (q) A G kl (q), (8) 

ki 

AG k (q) = G« +q/2 - Gt n /2 , AS k (q) = £« +q/2 - ££_ q/2 (9) 

where E k is the self energy and £7 kk /(q) is the irreducible 4- vertex appearing in the Bethe- 
Salpeter equation. This vertex has the diffusion pole for k + k' — > [1, 2] 

^Kq)=^ q ) + _^ +o fjy; k q) )(k+k0 , no) 

due to time-reversal invariance. The left-hand side of Eq. 8 is regular at the transition 
point, whereas the integrand on the right-hand side diverges as 1/uj in the localized phase 
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for all k, k' in the uo — > limit. This singularity should be cancelled after integration over 
k', which involves D(u>,h + k.') and imposes stringent requirements on the approximation 
used for calculation of the spatial dispersion of the diffusion coefficient. 

2. The first attempt to deal with these problems was made by Vollhardt and Wolfle 
[1]. They used approximation for t/ kk '(q) suggested by the weak localization theory [7] 
(corresponding to U^(q) = const, F(k, k', q) = const in Eq. 10) and solved approximately 
the Bethe-Salpeter equation. Their solution corresponds to the following simple estimate. 
The quantity U^'lq) plays the role of a 'transition probability' in the quantum kinetic 
equation, and one can use the analogue of the r-approximation, D oc I oc (U)' 1 (I is the 
mean free path, (...) denotes averaging over the momenta), to obtain the relation 

D~ const (u + F [ , ^ — ) , (11) 

coinciding with the self-consistency equation of the Vollhardt and Wolfle theory. With the 
increase of disorder, the 'transition probability' has anomalous growth due to diminishing 
of the diffusion coefficient and provides the possibility for vanishing of the latter. If the 
spatial dispersion of D(uj,q) is neglected, Eq. 11 becomes closed and allows to determine 
the critical exponents for the conductivity a and the localization length £, 

<J ~ T S , f ~ T~ V (12) 

(r is the distance to the transition). Setting D = const (u) ~ a in the metallic phase and 
D ~ (— ioo)^ 1 in the localized phase, one can obtain 



1 , d > 2 ; 



' l/(d-2) , 2<d<4 

(13) 

1/2 , d > 4 



Thus, the Vollhardt and Wolfle theory [1] provides the adequate qualitative description of 
the Anderson transition. However, due to neglecting the spatial dispersion of D(u,q) the 
problem of the order parameter symmetry was overlooked, while approximations U^,(q) = 
const, F(k, k', q) = const lead to the crude violation of the Ward identity. 

The problem of the spatial dispersion of the diffusion coefficient manifests itself also in 
the a-models formalism [8], which is usually considered as the most rigorous approach to 
the localization theory. These models are derived by using the saddle-point approximation 
for integration over 'hard' degrees of freedom in the functional integral and subsequent 
expansion in gradients. The 'minimal' a-model contains only the lowest (second) powers of 
gradients and corresponds to neglecting the spatial dispersion of the diffusion coefficient. 
One can suspect that it violates self-consistency of the theory. Attempt to take the spatial 
dispersion of D(u>,q) into account and include the terms with higher gradients into the 
Lagrangian of the cx-model leads to anomalous growth of these terms in the course of the 
renormalization-group transformations [9] : the so called 'high-gradient catastrophe' occurs. 
The stable fixed point, corresponding to the true spatial dispersion, is still not found. 
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The problem of the Ward identity was completely solved in the paper [2]. The crucial 
point consists in the use of the spectral representation for the quantum collision operator L 
(which is the quantum analogue of the Boltzman collision operator). Then cancellation of 
the 1/uj singularity in the right-hand side of Eq. 8 follows from the orthogonality relations 
for eigenvectors of L. At present, the spectral representation for the operator L looks as 
a necessary ingredient of any consistent theory of the Anderson transition: it is the only 
place where approximations can be safely made without violation of the Ward identity. 

The paper [2] deals also with the general problem of evaluation D(uj,q) near the An- 
derson transition. It appears, that the spectrum of the operator L possesses the non-trivial 
hierarhial structure. The condition of stability of this hierarhy in respect to infinitismal 
perturbations can be expressed in a form of a self-consistency equation, which replaces 
the crude Eq. 11. Solution of this equation is seeked for an arbitrary form of the spatial 
dispersion of D(u,q), but self-consistency is reached only for the solution with the weak 
g-dependence. As a consequence, estimation of the integral in Eq. 11 remains unchanged 
and and we return to the result (13) for the critical exponents. 

The paper [2] takes into account only evident symmetry of the system and considers 
a situation of the general position, compatible with this symmetry. Such reasoning is 
typical for the mean field theory and may be wrong due to existence of hidden symmetry. 
However, there are arguments that the theory [2] is something more than the mean field 
theory. Indeed, the existence of hidden symmetry is expected only for the critical point 
itself; as a result, the mean- field theory does not give the true critical behavior, but correctly 
describes the change of symmetry. From viewpoint of the change of symmetry, the Anderson 
transition is similar to the Curie point for an isotropic n-component ferromagnet in the 
limit of n — > oo: it is related with simultaneous vanishing of all coefficients in the expansion 
of D(0, q) over q 2 . The model of an infinite-component ferromagnet is exactly solvable [10] 
and its critical exponents appear to be in complete agreement with the results (13) of the 
straightforward analysis. This is an argument, that the symmetry of the critical point is 
established correctly and that the critical exponents are determined exactly. 

Another argument is given by the fact that Eq. 13 agrees with practically all reliable 
analytical results obtained for specific models [2, 11], and with experimental results s ~ 1, 
;/ k 1 for d = 3 [12, 13] obtained independently in the measurements of conductivity 
and the dielectric constant. It agrees also with the one-parameter scaling hypothesis [14], 
though not one word about scaling was said in derivation of the self-consistency equation. 
Unfortunately, the results (13) are in a serious contradiction with numerical experiments 
[15], though the latter have their own problems [16], being in conflict not only with (13) 
but also with some rigorous theorems. 

In conclusion, a scenario of the Anderson transition suggested in Ref. 2 cannot be con- 
sidered as rigorously established but has the good perspectives to be exact. At present, it 
is the only scheme which makes it possible to analyse the spatial dispersion of the diffusion 
coefficient near the critical point. 

3. In the current literature, behavior of the diffusion coefficient near the Anderson 
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transition is discussed in terms of one-parameter scaling [17], and the following ansatz is 
used for a critical point [18] 

D(uj,q)~ U j 1, / d q d - 2 - T >. (14) 

The exponents in Eq. 14 are chosen in such way that the static diffusion coefficient D L of 
a finite system of size L, determined from D(u>, q) by relation 

D L ~D(D L /L 2 ,1/L), 

has a behavior D L ~ L 2 ~ d predicted by the one-parameter scaling theory [14]. Relation 
i] = d — D 2 was claimed [19], connecting the exponent rj with the fractal dimension D 2 of 
wave functions. Numerical results for r\ and D 2 in three dimensions have a large scattering 
( D 2 = 1.7±0.3 [20], D 2 = 1.6±0.1 [21], D 2 = 1.33±0.02 [22], r] = 1.2±0.15, r] = 1.3±0.2, 
rj = 1.5 ± 0.3 [19], D 2 = 1.68 [23], D 2 = 1.30 ± 0.05 [24], D 2 = 1.28 ± 0.02 [25]) but in 
general show violation of an equality rj — 1, which follows from the analysis of Ref. 2. It is 
possible to consider this fact as an evidence of hidden symmetry neglected in [2]. However, 
few critical remarks can be made in relation with the cited researches. 

(a) The analysis in [18, 19] is based on the relation between the Fourier transform of 
the density correlator (5) and the diffusion coefficient 

u 2 + [D{u, q)q 2 \ 2 

which is valid only for real D(cu,q). The general relation can be found in [2] and is given 
there by Eq. 31: 

S(u,q)~hm D{U ^ q \ 2 . (16) 
V ,H! lo + D{u,q)q 2 V ; 

It is clear from the above discussion that the low frequency diffusion coefficient is real in the 
metallic state and pure imaginary in the localized phase, while a complicated rearrangement 
of its analytical structure occurs in the vicinity of the transition point. Thus, the most 
delicate feature of the Anderson transition is completely ignored in [18, 19]. 

In fact, only one real quantity S(u,q) is measured in numerical experiments, and one 
cannot extract two functions R,eD(u,q) and Im D(u,q) from such a measurement 1 . We 
see, that numerical experiments do not provide any direct information on the diffusion 
coefficient. All conclusions concerning D(u,q) are based on strong assumptions and cannot 
be considered as reliable. 

(b) In our opinion, there is some problem with the relation 77 = d — D 2 . According 
to Wegner [26], D 2 = 2 — e in the (2 + e)-dimensional case, suggesting rj = 2e. However, 
Wegner's result was obtained for the 'minimal' a-model, where the spatial dispersion of 
D(u,q) is neglected (see Sec. 2 above); it looks as internal inconsistency. 



1 In principle, one can use the Kramers-Kronig relation and obtain a closed system of equations for 
these quantities, but the contemporary algorithms are far from fulfilling this program. 
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(c) The fractal dimension D 2 is determined by the behavior of the participation ratio 
P as a function of the system size L, P ~ L D2 . Theoretically, such behavior is expected 
only for the critical disorder W c , but in practice the power law dependences P ~ L a are 
observed for wide range of disorder [21] (see a detailed discussion in [27]). In order the 
measurement of D 2 was possible, the value W c is taken from other experiments, which are 
not directly related with the properties of wave functions. Analysis of the contemporary 
situation shows [16] that essential revision of the conventional results for W c is practically 
enevitable. As a result, a value D 2 can change essentially. 

4. Relation of the diffusion poles with the Anderson transition was discussed recently 
in the series of papers by Janis and Kolorenc [28-31], where the following conclusions were 
reached: 

(i) Existence of the diffusion pole in the localized phase is incompatible with the Ward 
identity; 

(ii) The diffusion pole is absent in the localized phase; 

(iii) The Ward identity is violated due to the averaging procedure. 

The initial point for these conclusions is the solution for the two-particle Green function 
obtained for high dimensionalities of space d, which (according to the authors) is asymp- 
totically exact in the parameter d^ 1 [28]. This solution does not posess the diffusion pole 
in the localized phase and does not satisfy the Ward identity. Instead to seek an origin 
of this surely defective result, the authors of [28-31] become anxious to defend it. Using 
reasonable, but unjustified form of {7kk'(q), they came to conclusion on the principal im- 
possibility to satisfy the Ward identity in the localized phase [29]. Its violation is argued 
by possibility that the Hilbert space of the eigenfunctions becomes incomplete due to the 
averaging procedure [30]. 

We have specially stressed in Sec. 1, that in order to satisfy the Ward identity one 
needs very delicate approximations: they should be not only asymptotically exact in some 
parameter but should be made 'in the proper place' 2 . We see nothing unusual that the 
results for high dimensions are not self-consistent, even if they are really correct in the 
leading order in d^ 1 . Indeed, the formal expansion of the right-hand side of Eq. 8 in powers 
of <i _1 may contain the well defined zero order term, while divergent integrals may occur 
in the higher orders. Consequently, the result of [28] for high dimensions does not provide 
the reliable foundation for further conclusions. 

The statement (i) on the principal impossibility of satisfying the Ward identity in the lo- 
calized phase (based on approximate structure of C/kk'(q)) contradicts the rigorous analysis 
of the paper [2] . The 1/ou singularity in the right-hand side of Eq. 8 is enevitably cancelled 
if the spectral properties of the quantum collision operator L are taken into account. 

2 There is nothing new in these considerations. It is well known that the first correction to the self energy 
has essentially more sense that the first correction to the Green function; the one-loop approximation for 
the Gell-Mann - Low function is more effective than the one-loop approximation to the invariant charge, 
and so on. 
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The statement (ii) on the absense of the diffusion pole in the localized phase contradicts 
the existence of the Berezinskii-Gorkov singularity which can be obtained not only 
in the framework of the self-consistent theory [1] (which is cited in [28-31]) but by direct 
analysis of the density correlator [3, 4, 5, 2] and in the framework of the instanton approach 
[4, 32]. This singularity is closely related with the Anderson criterion of localization [33] 
in modification by Economou and Cohen [34] (see [4] and [32]), which can be rigorously 
proved for ID systems [35]. In fact, the existence of the 1/u singularity can be easily seen 
for a zero- dimensional case, which corresponds to the Anderson model on a single site. The 
exact Green function 

after trivial averaging takes a form 

— oo 

(P(V) is a distribution of V), giving a density of states 

u{E) = --Im {G R (E)} = P(E) . (19) 

7T 

Then the quantity (p RA , 

oo 

0- = < G *(£ + «)G* ( E)) = / E _ v+ .„ ^317-5 P(V) "V , (20) 



— OO 



has the 1/uj singularity for uj — > 0, 



0- = ^>, (21) 

— 



in accordance with Eq. 3 (the momentum dependence is absent in the zero-dimensional 
case, and all momenta are supposed to be equal to zero; so q = and A(0) = 1). The 
results (17-21) are also valid for a space of arbitrary dimensionality in the limit of strong 
disorder. Already these simple results are not reproduced by sophisticated expressions in 
[28-31]. 

It is well known, that the coherent potential approximation (CPA), considering as for- 
mally exact in the leading order in has an essential qualitative defect [36]: it does 
not describe the fluctuational tail of the density of states, which is exponentially small 
and cannot be found in any finite order of perturbation theory. This tail can be explicitly 
calculated by the instanton method both for d < 4 [32, 37, 38] and for higher dimensions 
[39, 11] 3 . According to [4, 32], the Berezinskii-Gorkov singularity has also nonperturbative 
origin, and it is naturally not reproduced by the extention of CPA used in [28-31]. 



3 Absence of instantons for d > 4 [40] should not be considered as the absence of localized states [41]. 
The theory is not renormalizable for d > 4 and should be considered on the lattice. The lattice instantons 
always exist. 
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The statement (iii) on the possible violation of the Ward identity due to the averaging 
procedure [30] is based mainly on philosophical considerations. In principle, such argu- 
mentation may be correct in some situations, but it has nothing to do with the present 
case: the Ward identity (8) was diagrammatically derived in [1] directly for the averaged 
quantities. Of course, one can discuss the possible nonperturbative contributions, but at 
the present case there is no ground for them: validity of the Ward identity corresponds to 
the physical requirements, and its violation is not expected. In fact, the mystical fear of 
nonperturbative contributions is not grounded: such contributions can be safely extracted 
from the diagrammatic expansions [11] if a correct procedure for summing of divergent 
series is chosen [42]. 

One can see, that all statements (i),(ii),(iii) are incorrect. The authors of [28-31] were 
wrecked on the 'hidden rocks', which were discovered many years ago. 

This work is partially supported by RFBR (grant 03-02-17519). 
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